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We study the nonlinear energy transfer around the peak of the spectrum of surface gravity waves 
by taking into account nonhomogeneous effects. In the narrow-banded approximation the kinetic 
equation resulting from a nonhomogeneous wave field is a Vlasov-Poisson type equation which in- 
cludes at the same time the random version of the Benjamin-Feir instability and the Landau damping 
phenomenon. We analytically derive the values of the Phillips' constant a and the enhancement 
factor 7 for which the narrow-banded approximation of the JONSWAP spectrum is unstable. By 
performing numerical simulations of the nonlinear Schrodinger equation we check the validity of the 
prediction of the related kinetic equation. We find that the effect of Landau damping is to suppress 
the formation of coherent structures. The problem of predicting freak waves is briefly discussed. 
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In many different fields of nonlinear physics, local 
nonlinear effects such as the modulational instability 
(MI) have played a very important role in spectral en- 
ergy transfer processes (see for example Jl-§1). For 
ocean gravity waves, the same instability (now commonly 
known as the Benjamin- Fcir instability) has been discov- 
ered independently by Benjamin and Feir [Q and by Za- 
kharov Q in the sixties. The instability predicts that 
in deep water a monochromatic wave is unstable under 
suitable small perturbations. In this framework it has 
been established that the MI can be responsible for the 
formation of freak waves ■ The aim of this Letter is 
to discuss a statistical approach that allows one to pre- 
dict the appearance of the modulational instability in a 
continuous spectrum. The theory is applied directly to 
ocean waves but it could as well be applied to many other 
nonlinear wave systems. 

In order to approach statistically the nonlinear energy 
transfer processes involved in the MI, one is interested 
in finding a suitable kinetic description. As far as ocean 
waves are concerned, today the most common models 
for wave forecasting (WAM, SWAM and newer genera- 
tion models, see Q) are based on the nonlinear energy 
transfer process that is ruled by the kinetic wave equation 
that has been derived independently by K. Hasselmann 
jl0| and by V. Zakharov ||. The theory predicts that 
the energy is transferred in an irreversible manner un- 
der the resonant 4-wave interaction. Besides the quasi- 
gaussian approximation |fllf , one of the major hypothesis 
required to derive this kinetic equation is the homogene- 
ity of the surface, i.e., < A(k)A*{k') >= n(k)5(k - k'), 
where A is a complex wave amplitude describing the evo- 
lution of the wave train, k and k' are wave numbers and 
n(fc) is the spectral density function. Unfortunately, the 
Hasselmann-Zakharov kinetic theory is not able to pre- 
dict the Benjamin-Feir instability. This basically means 
that the kinetic equation is unable to predict local nonlin- 



ear effects such as freak waves. Indeed the MI is the result 
of an interaction of waves that are phase-locked: the car- 
rier wave is phase-locked with the side-bands and there- 
fore the mechanism cannot be predicted by the assump- 
tion that the Fourier components are delta-correlated. 
Thus, if the hypothesis of homogeneity of the system 
is relaxed, an improved kinetic equation can be derived 
which is able to show a random version of the Benjamin- 
Feir instability. Actually, for surface gravity waves, this 
improvement is contained in the pioneering wor k by Al- 
ber |fl2|] , followed by the works of Crawford et al. j[3[ and 
Jansse n |1l4| ,p"5|| . The prediction of the theory developed 
in [fl2| - |l4| has never been verified numerically or experi- 
mentally. Unfortunately these works, written more than 
15 years ago, did not receive the deserved attention by 
the wave forecast /hindcast community and the ideas de- 
veloped remained basically unconsidered with the only 
exception of the recent work by M. Stiassnie pq| . More 
recently, but independently of the works |12|-|l4|^ a simi- 
lar approach has been developed for the large-amplitude 
electromagnetic wave-envelope propagation in nonlinear 
media |l7f , for the quantum- like description of the longi- 
tudinal charged-particle beam dynamics in high-energy 
accelerating machines in the presence of an arbitrary 
coupling impedance [jl8| and for the resonant interaction 
between an instantaneously-produced disturbance and a 
partially incoherent Langmuir wave pl| . 

In this paper, we outline the approach formulated in 
p2| |L4| and discuss the modulational instability for ran- 
dom wave spectra. In particular we identify the values 
of the parameters of the JONSWAP spectrum for which 
the spectrum itself is unstable. Moreover, using numer- 
ical simulations we associate the MI with the presence 
of coherent structures in physical space, pointing out the 
importance of the role played by the phenomenon of Lan- 
dau damping. 

We start by considering the simplest and prototypical 
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weakly nonlinear equation for water waves: the Nonlin- 
ear Schroedinger equation (NLSE). This choice is dic- 
tated by the following motivations. First of all, we want 
to avoid dealing with the dynamics of the four-wave res- 
onant interactions because otherwise we would not be 
able to discern whether the effects we are describing are 
due to inhomogeneity or are already contained in the 
Hasselmann-Zakharov theory. We point out that, while 
the two dimensional NLSE (and higher order equations 
(see 20 )) contains the four- wave resonant interaction 
mechanism plfl , the one-dimensional NLSE does not in- 
clude this sort of energy transfer, because of the shape 
of the linear dispersion relation (see for instance fp2[). 
Second, while the NLSE should not be appropriate for 
describing the dynamics in the tail of the spectrum or 
in the inertial range, it should describe with satisfactory 
accuracy the short time scale behavior around the peak. 
Formally the NLSE can be derived from the Zakharov 
equation [^| under the narrow band approximation; in 
dimensional form the equation reads: 



OA 
~dt 



d 2 A 

dx 2 



-w\A\ 2 A = 0, 



(1) 



where in deep water /i — cjo/8fc§ and v = took 2 ,/ 2, with 
too the carrier angular frequency and fco the respective 
wave number. 

Eq. ([!]) is our starting point for deriving the required 
kinetic equation. Following Alber the Wigner-Moyal 
transform |23| can be applied directly to the NLSE. This 
transform allows one to give a representation of a func- 
tion A(x) both in configuration space, x, and in phase 
space or wave number, fc, namely 



n(x, k) 



2tt 



< A*(x + y/2)A(x - y/2) > e~ tky dy 



(2) 



Alternatively, one may follow the approach presented 
both in and in consisting in writing the NLSE 
in Fourier space and then writing an evolution equation 
for the correlator n(k,k') =< A(k)A*(k') >. Both ap- 
proaches need a closure in order to reduce the fourth- 
order correlator into the sum of the product of the two- 
point correlation functions. Using any of the methods 
outlined the resulting kinetic equation is the following 
von Neumann- Weyl-like equation 

dn(x,k,t) , 2fik dn(x,k,t) { ^ 



dt 



dx 



(2m+ l)!2 2m+1 



771=0 



d 2m+i < \A(x,t)\ 2 > d 2m+1 n(x,k,t) 



dk 2m+1 



with 



< \A(x,t)\ 2 >= / n(x,k,t)dk. 



(3) 



(4) 



If the limit for small k is taken the equation resembles the 
Vlasov-Poisson equation in plasma physics that is well 



known to describe the Landau damping phenomenon. 
The standard way to proceed consists in letting the dis- 
tribution function n(x,k,t) and the field A(x,t) be ex- 
pressed in terms of an equilibrium value plus a small 
perturbation and study the dispersion relation of the lin- 
earized equation for the perturbation. After standard 
algebra the following dispersion relation is obtained (see 
also fEMl): 



v f n (k + K/2)-n (k-K/2) 

where no is the homogeneous envelope spectrum. 

According to experimental studies, the surface wave 
spectrum for gravity waves is given by the JONSWAP 
spectrum: 



P(k) = 



— p 2 L fc J 



2fc 3 



7 



(vTfc)-y(fco)) 2 

25' 2 fc 



(6) 



with a, 7 and 5 constants (S is usually set to 0.07, while 
a and 7 depend on the state of the ocean). In order 
to be consistent with the NLSE, as is done in JjUfl , we 
Taylor-expand the spectrum around its peak and obtain 
the following Lorcnzian spectrum: 



P(k) 



H 2 



P 



IQirp 2 + (fc - fc ) 2 



where 



P = 



8k 2 S 2 



24S 2 + Login) 



and H„ = 4 



~2E 3 / 2 kl 



(7) 



(8) 



H s is the significant wave height calculated as four times 
the standard deviation of the wave field and p corre- 
sponds exactly to the half-width at half-maximum of the 
spectrum. It can be shown ||l3|] that for a symmetric 
spectrum P(k) of the surface elevation, the spectrum for 
the complex envelope is rto(fc) = 4P(fc — ko), therefore a 
factor of four must be taken into account. Substituting 
in (^J) we obtain the following dispersion relation: 



Q = K{y/K 2 n 2 - H 2 v^ - 2inp). 



(9) 



If K 2 < H 2 v/fi, the first term on the right hand side is 
responsible for the MI (note that in the limit as p — > 0, 
the dispersion relation <M) gives the Benjamin-Feir insta- 
bility). The last term on the right hand side is responsi- 
ble for the Landau damping phenomenon p5| . Therefore 
there is a competition between exponential growth and 
damping of the perturbation that depends on the param- 
eters a and 7 of the Lorenzian (or JONSWAP) spectrum. 
If p > y/ — K 2 ii/A + H 2 v/A the damping dominates the 
MI, the opposite will occur if p < y/—K 2 /i/A + H^v/A. 
In Fig. [l] we show the instability diagram in the a — 7 
plane. In general spectra with higher values of a and 7 
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are more likely to show the MI. This results gives sup- 
port to the results in (§1 where it has been found that 
higher values of a and 7 increases the probability for the 
formation of freak waves. 
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FIG. 1. Instability diagram in the a — 7 plane. 

Before comparing numerical simulations with the pre- 
diction of the Wigner-Moyal-like kinetic equation, we 
would like to point out that our analysis has been carried 
out linearizing around an equilibrium state of the system. 
In fact, the dispersion relation (|^) has been derived for 
small amplitude perturbations. However, in the natural 
long-time evolution of a nonlinear wave train, perturba- 
tions are in general not small. Consequently, perturba- 
tions in the simulations should evolve according to the 
governing equations. To this end, numerical experiments 
have been carried out allowing the nonhomogeneities to 
develop by themselves according to the nonlinear dynam- 
ics of the NLSE. 

Numerical simulations of Eq. (|l|) have been computed 
using a standard pseudo-spectral Fourier method. Initial 
conditions for the free surface elevation £(cc, 0) have been 
constructed as a random process pr| , i.e. a linear super- 
position of Fourier modes (|^) with random phases. The 
Hilbert transform is used in order to convert the free sur- 
face £ to the complex envelope variable A of the NLSE. 
The dominant wave number for the numerical simulation 
was selected to be fco = 0.1 m . This last choice is not 
restrictive: the parameters that rule the dynamics in the 
spectrum are its width and the steepness which, once fco 
is fixed, are univocally determined by a and 7. 

We start the discussion on the numerical results by 
showing the evolution of in the x — t plane for 

a case in which 7 = 3 and a = 0.02 (p/fco — 0.25 and 
H s k /2 ~ 0.14 ), Fig. |. According to Fig. [I] the spec- 
trum should be unstable. How is this instability man- 
ifested in configuration space? From Fig. ^ we note 
the presence of a "coherent structure" , i.e. a structure 
(oblique darker zones in the x — t plane) that persists in 
the presence of nonlinear interactions and maintains sta- 



tistically its shape and velocity during propagation (note 
that periodic boundary conditions are used). The pres- 
ence of such coherent structures is related to the values 
of a and 7 and not to the particular random phases se- 
lected. 
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FIG. 2. |A(x,t)| from numerical simulation of the NLSE. 
The initial condition is characterized by a Lorentian spectrum 
with 7 = 3, a — 0.02. A coherent structure is evident in the 
x — t plane. 
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FIG. 3. from numerical simulation of the NLSE. 

The initial condition is characterized by a Lorentian spectrum 
with 7 = 1, a = 0.02. The field in the x — t plane appear to 
be random with out any evidence of coherent structure. 

Every random realization shows similar results even 
though the resulting coherent structures may have differ- 
ent velocity and amplitude. The nonlinear stage of MI is 
therefore responsible for the formation of such coherent 
structures in the x — t plane. Indeed it is possible to show 
that the NLSE has periodic solutions such as, for exam- 
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pie, breathers or unstable modes [|27j,|28|. These solutions, 
which are linearly unstable, are nevertheless very robust. 
Moreover they can grow up to more than three times 
the initial unperturbed solution and therefore have also 
been addressed as simple models for freak waves [^8],^9| . 
In contrast to solitons that have constant amplitude in 
time, these unstable modes are characterized by a contin- 
uous exchange of energy among the Fourier modes. The 
energy is transferred from one mode to another and back 
again: the process is completely reversible and therefore 
coherent structures persist in physical space. We stress 
that these kinds of solutions appear naturally from ini- 
tial conditions with random phases. We now consider 
a case with 7 = 1 and a = 0.02 (p/ko ~ 0.58 and 
H s ko/2 ~ 0.12); note that decreasing 7 has the effect 
of making the spectrum broader. This case, according 
to the stability criteria (see Fig. [j]), should not present 
any instability because of the effect of Landau damping. 
We show the output of our numerical simulation in Fig 
|3|. In the x — t plane, as a confirmation of what we have 
previously stated, there is no evidence of any coherent 
structure; numerical simulations with different random 
phases are in accordance with the result just shown. 

In conclusion, the numerical simulations show that the 
presence of the coherent structures in the x — t plane is 
related to the instability of the wave field. The Landau 
damping phenomenon suppress the instability and pre- 
vent the formation of coherent structures. Many physical 
questions remain open. For example it would be inter- 
esting to investigate the case of a two dimensional wave 
field. It is well known that the NLSE in 2+1 is not inte- 
grable and the dynamics of coherent structures is still far 
from being understood. Numerical simulations with the 
fully nonlinear Euler equations are also under consider- 
ation in order to extend the validity of the result. Con- 
cerning the problem of wave forecasting, we may state 
that if one is interested in predicting freak waves arising 
from the MI, this new form of the kinetic equation, which 
includes nonhomogeneous effects should be considered. 
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